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Abstract 

The potential between two heavy-light mesons as a function of the heavy 
quark separation is calculated in quenched SU(3) lattice QCD. We study 
the case of heavy-light mesons with a static heavy quark and light quarks of 
mass close to the strange quark mass. We explore the case of light quarks 
with the same and with different flavours, classified according to the light 
quark isospin. We evaluate the appropriate light quark exchange 
contributions and explore the spin-dependence of the interaction. 
Comparison is made with meson exchange. 

PACS numbers: 13.75.Lb, ll.15.Ha, 12.38.Gc, 25.80.-e 



1 Introduction 

The progress in lattice QCD has so far been mainly restricted to systems of 
three quarks or less. However, there is also considerable interest in obtaining 
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predictions from first principles for multi-quark systems which can be decom- 
posed into more than one colour singlet. In addition to the complicated cases 
of nuclei, simple multi-quark systems have been proposed to exist as bound 
states [1, 2, 3]. Four quarks forming colour singlets or as bound states of two 
mesons are candidates for particles lying close to meson-antimeson thresh- 
old, such as ao(980), /o(980) {KK), /o(1500), /2(1500) {uju, pp), /j(1710) 
{K*K*), ^(4040) {D*D*), T(10580) {B*B*) [4]. 

Systems with heavy quarks should be more easily bound provided the poten- 
tial is attractive, since the repulsive kinetic energy of the quarks is smaller, 
while the attractive two-body potential remains the same. In so-called deuson 
models [5] the long-range potential between two mesons comes from onc-pion 
exchange, suggesting that meson-meson systems are significantly less bound 
than meson-antimeson systems. Other models used for four-quark systems 
include string- flip potential models (see Ref . [6] for a review) , bag models [2] , 
and a model- independent approach [7]. Four-quark states with two heavy 
quarks have been predicted to be stable [8]. Most models give stability for 
systems where the heavy quarks have the b mass, but long range forces might 
push the required heavy-to-light mass ratio down so that ccqq states would 
be bound as well. 

Static four-quark systems [9] (and references therein) have been previously 
studied for a set of geometries representative of the general case and a model 
was constructed that reproduces one hundred ground and excited state ener- 
gies with four independent parameters [10]. The model is based on ground- 
and excited state two-body potentials and multi-quark interaction terms. 
The results show that, in a two-body potential approach to understanding 
multi-quark interaction, the effect from gluonic excitations is needed, and 
their relative contribution to the binding becomes more important (even 
dominant) at larger distances. Flux distributions corresponding to the bind- 
ing energies of four static quarks are studied in Ref. [11]. 

Moving on to more realistic systems, we now study in detail the potential 
between two heavy-light mesons. Exploratory studies of two-meson systems 
have been made for the cross diagram only (Fig. 1 below) for SU(3) colour [12] 
and for both diagrams in Refs. [13], [14] for SU(2), SU(3) colour respectively. 

We take the mass of one quark in each meson to be heavy - the prototype 
being the B meson. This is in the spirit of the heavy quark effective theory 
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approach which describes the leading term (the static hmit) and the correc- 
tions of higher orders in l/mq. In the static approximation for the heavy 
quarks, the pseudoscalar B meson and the vector B* meson will be degener- 
ate - whereas they are split by 46 MeV experimentally. Since we shall often 
have occasion to treat this degenerate set together we describe this 
the B meson. In analogy to the Born-Oppenheimer approximation, we will 
then discuss the potential energy between static B mesons. 

For the light quarks, we use the full relativistic description with a fermion ac- 
tion which is the O(a^) improved Sheikholeslami-Wohlert clover action with a 
tadpole-improved coefficient. We should in principle evaluate the interaction 
for several light quarks masses and then extrapolate to the physical values. 
In this preliminary study, we fix the light quark mass at around the strange 
mass. We do however consider the case of two flavours of quark - so allowing 
a discussion of different isospin states. The main reason why this study is 
difficult to perform on a lattice is that the light quark propagators are needed 
from many different sources. To achieve this we make use of the technique 
of evaluating the light quark propagators as stochastic estimates [15] using 
maximal variance reduction introduced in Ref. [16]. 

Quenched lattices are used with SU(3) colour and static heavy quarks with 
light quarks of approximately the strange quark mass. Preliminary versions 
of this work have appeared [21]. Here the isospin and spin degrees of freedom 
are discussed in detail. We compare our results for small separation R with 
the known spectrum of baryons with one heavy quark ( A5 and S5) . This will 
enable us to discover if a heavy diquark is a good description. Note that 
this link which we find to baryons at small separation R cannot be explored 
using SU(2) of colour. We also compare our results with the expectations of 
meson exchange. We find that at larger R, this is a useful guide to the inter- 
action strength and, for pion exchange, we are able to make a quantitative 
comparison. We comment on the agreement with other models, one of them 
being the potential model for static systems applied in this more dynamic 
case [19]. 
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2 BB interactions in the static approximation 



We take the mass of one quark in each meson to be very heavy - the prototype 
being the B meson. The static hmit is then the leading term in the heavy 
quark effective theory for a heavy quark of zero velocity and there will be 
corrections of higher orders in 1/ rriQ where mq is the heavy quark mass. In 
the limit of a static heavy quark, the heavy quark spin is uncoupled since 
the relevant magnetic moment vanishes which implies that the pseudoscalar 
B meson and the vector B* meson will be degenerate. This is a reasonable 
approximation since they are split by 46 MeV experimentally, which is less 
than 1% of the mass of the mesons. Since we shall often have occasion to 
treat these two mesonic states as if they were degenerate, we describe them 
collectively as the B meson. Because of the insensitivity to the heavy quark 
spin, it is then appropriate to classify these degenerate B meson states by 
the light quark spin: so there are only two independent spin states. The 
system of two heavy-light mesons at spatial separation R will be referred to 
as the BB system. With both heavy-light mesons static, this BB system is 
described by the spin states of the two hght quarks in the two mesons. Thus 
there are four possible states and we need to classify the interaction in terms 
of these spin states. 

This situation is very similar to that of the hydrogen molecule in the Born- 
Oppenheimer approximation - with, however, the additional possibility that 
the two 'electrons' can have different properties. Another similarity is with 
the potential between quarks which has a central component and then scalar 
and tensor spin-dependent contributions. 

Each B meson will have a light quark flavour assignment. For the BB system, 
it will be appropriate to classify these states according to their symmetry 
under interchange of the light quark flavours. For identical flavours (eg. 
ss or uu), wc have symmetry under interchange, whereas for non-identical 
flavours (eg. su or du), we may have either symmetry or antisymmetry. For 
two light quarks, it is convenient to classify the states according to isospin 
as / = 1 (with uu, ud + du and dd) or I — (with ud — du). 

We now present a discussion of the possible states of two B mesons. As a 
guide we show in Table 1 the states for the case of an S-wave BB system in 
the limit of static B mesons. Wc must have overall symmetry of the wave- 
function under interchange and, assuming symmetry for spatial interchange. 
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the flavour, total light quark spin (Sq) and total heavy quark spin (Sh) must 
be combined to achieve this. Thus in the limit of an isotropic spatial wave- 
function, there will be the four different ground state levels of the BB system 
as shown in Table 1 since the three states with different but the same 
light quark isospin Ig and spin Sg will be degenerate in the static limit. We 
will label these states by Jg, Sg for subsequent discussion. We also show 
which physical B and B* mesons couple to these states. This table can also 
be extended to L 7^ levels. In particular, we shall later see that a tensor 
interaction may be present, in which case the Sq — 1 ground states will show 
an admixture of L = and of L = 2. 

Table 1: Allowed BB states with L = 

Ig Sg Sb BB BB* B*B* 

1 i i 0+ Yes Yes 
1111+ Yes 
1112+ Yes 

1 0+ Yes Yes 
10 1+ Yes Yes 

11+ Yes Yes 

When R = 0, the situation is special since the colour of the two static quarks 
can be combined. This net colour can be in a anti-triplet (antisymmetric 
under particle exchange) or a sextet (symmetric under particle exchange). 
The former case is just that of the static baryons. This equivalence implies 
that the Iq — 1, Sq — 1 state will have the same light quark structure as 
the Eft baryon, while the Iq — 0, Sq = state will be as the baryon. The 
other two allowed BB states a,t R — correspond to a static sextet source. 

In a Born-Oppenheimer treatment of the BB system, we will need to consider 
the potential energy for the B mesons at rest at separation R. This BB system 
can be classified under rotations about the separation axis, here taken as the z 
axis, and under interchange of the two mesons. Taking the z axis to quantise 

the light quark spin, we have states with — ±1, namely | + +), | ) 

and with = 0, namely | H — ) ± | — h). Since is conserved, we can 
discuss the interaction energy in terms of a triplet state Et corresponding 
to the J2 = ±1 cases, and then the = sector can be described by a 
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singlet state with energy Eg between initial and final states (| H — ) — | — h)) 
and by another triplet state for initial and final states (| H — ) + | — h)). 
These three energies can be related to a more conventional treatment using a 
central, spin- dependent and tensor potential. We shall instead mainly focus 
on the singet, Es, and triplet averaged over orientations, {2Et + E'j) /?>. Since 
the heavy quark spin docs not interact, the symmetric and antisymmetric 
combinations constructed from the heavy quark spin will allow any overall 
symmetry under interchange for the overall spin assignment. 

In our actual numerical calculation, we use a relativistic treatment of the light 
quark spin but in the context of a static heavy quark with Dirac propagator 
structure (1 + 74). This enables us to simplify the Dirac 7-matrix algebra 
between initial and final B mesons (created by ^756) and B* mesons (created 
by (pfih). This approach also leads to three independent observables which 
we determine as 

= ((++) + (—))■((++) + (— )) 

<^.(^) -((++)-(—))•((++)-(— )) 

Cs{x) = Csiy) = (+-).(-+) + (-+).(+-) 

with notation (13). (24) for -B1-B2 — ^ -B3-B4 with the sign of the light quark 
spin (5*12, etc) given. 

In practice the observable given above by Cs{z) is also evaluated with the 
spatial separation Rinx and y directions which gives an equivalent method to 
obtain Cs{x) and Cs{y). By symmetry, the latter two observables are equal 
on average. Note that the BB — > BB correlation is given by C/, whereas 

BB* ^ B*B is given by C,. 

It is not sufficient just to look at processes such as BB BB since, in the 
heavy quark limit, there will also be other channels such as BB — > B*B* which 
are coupled. We then analyse the matrix of correlations between all such 
channels and find the basis that diagonahses it. This leads to certain linear 
combinations of correlations which describe these diagonal elements. This 
explicit fcrmionic approach must reprodiicc the conclusions reached above 
by using the heavy quark limit. The relationship between these approaches 
is that, at large t, 

Cj + Cs{z) ^ e-^-* (1) 
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Ci - Cs{z) - Cs{x) - Cs{y) -> e-^^* (2) 

Ci - Cs{z) + Cs{x) + Cs{y) ^ e-^T* (3) 

It turns out that the same combinations (those given in the above equations) 
occur for both the case of symmetry under exchange of initial particles and 
for the case of antisymmetry. This can be understood in the heavy quark 
effective theory, as discussed above, since the combinations are in terms of 
the light quarks spins, leaving the heavy quark spins to be combined either 
in symmetric or antisymmetric states. 

The structure of the correlations to be evaluated, in terms of the light quark 
propagator G and the gauge product for the static line in the negative-going 
t-direction of U is then 

Ci{t) = (G^r(0, 0; 0, t)U''\Q) Gtlie^R, 0; e,R, t)W\e,R)) (4) 

where the colour indices a, 6, c, d and the Dirac indices j, /c, / are associated 
with vertices as in Fig. 1. The sum over Dirac indices is only from 1 to 2 
since the heavy quark has a spin projection factor. These contributions can 
be evaluated for every choice of origin on a lattice which is translationally 
invariant. For the spin-dependent part with component p, we have 

Cl{t) = {GfM 0; 0, Gf^aUe.R, 0; e,R, t)U'^{e,R)) (5) 

For the 'cross' diagram the colour and spin sums are different, for example 
the contribution to C/ is given by 

Ciit) = -{Gt^iO,0;e,R,t)U^''ie,R) G^^(e,i?, 0; 0, t)f/»^(0)) (6) 

where the negative sign comes from the Grassmannian nature of the fermions. 
For states symmetric under light quark interchange (eg. 7=1), then the 
sum of uncrossed and crossed diagram is needed, where the above minus 
sign is incorporated into the crossed diagram - this plays the role of the 
Pauli principle. For states antisymmetric under light quark interchange (eg. 
7 = 0), the difference of uncrossed and crossed diagrams is needed. 

3 Fermion formalism 

The diagrams we need to evaluate are illustrated in Fig. 1. We need light 
quark propagators from more that one source - so the conventional approach 
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a, i b, j 




Cy k dy I 

Figure 1: Diagrams showing the interaction between two B mesons: the hght 
quarks are shown as wiggly hues. 

of inverting from a single source is impractical. One feasible way forward is 
to use a stochastic inversion method which allows the evaluation of quark 
propagators from any site to any other site. The stochastic method has 
already been shown to be more efficient than the conventional inversion from 
one source for mesons made of heavy-light quarks [16], and it does allow the 
flexibility to evaluate the required combinations of correlations readily. For 
this reason it allows a thorough study of this area. 

Stochastic propagators [16, 15] are one technique to invert the fermionic ma- 
trix for the light quarks. They can be used in place of light quark propagators 
calculated with the usual deterministic algorithm. The stochastic inversion 
is based on the relation: 

G^J = M-l = ^Jl^<P {MJk<Pkr<P^ exp [-<p*{M^M)^j<Pj) (7) 

where, in our case, Al is the improved Wilson-Dirac fermionic operator and 

the indices i,j,k represent simultaneously the space-time coordinates, the 
spinor and colour indices. For every gauge configuration, an ensemble of 
independent fields (pi (we use 24 following [16]) is generated with gaussian 
probability: 

P[0] = lcxp(-0*(A^t^),,0,) (8) 
All light propagators are computed as averages over the pseudo-fermionic 
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samples: 

G 



'13 



or (9) 
, 75(0*(-M0)i)75 

where the two expressions are related by Gij — ^^G^jff^. Moreover, the maxi- 
mal variance reduction method is applied in order to minimise the statistical 
noise [16]. The maximal variance reduction method involves dividing the lat- 
tice into two boxes (0 < t < T/2 and T/2 <t <T) and solving the equation 
of motion numerically within each box, keeping the pseudo-fermion field on 
the boundary fixed. According to the maximal reduction method, the fields 
which enter the correlation functions must be either the original fields or 
solutions of the equation of motion in disconnected regions. The stochastic 
propagator is therefore defined from each point in one box to every point in 
the other box or on the boundary. For more than one propagator from one 
box to the other, we need to use different stochastic samples for each. This 
is completely analogous to the technique used to discuss the A;, meson [16]. 

The numerical analysis used 24 stochastic samples on each of 60 quenched 
gauge configurations, generated [16] on a 12^ x 24 lattice aX (3 = 5.7, corre- 
sponding to a^^ = 1.10 GeV. With improved clover coefficient Csw = 1-57, 
we use a value of Ki = 0.14077 which corresponds to a bare mass of the 
light quark around the strange mass and gives a pseudoscalar to vector mass 
ratio of 0.650(7). The chiral hmit corresponds to Kc — 0.14351 [18]. Error 
estimates come from bootstrap over the gauge configurations. We also used 
20 quenched gauge configurations on a 16'^ x 24 lattice to check finite size 
effects with the same parameters as above. Allowing for the self-averaging 
effect of the larger spatial volume, this data set has similar weight to that at 
the smaller volume. 

In smearing the hadronic interpolating operators, spatial fuzzed links are 
used. Following the prescription in [16, 17], the fuzzed links are defined 
iteratively as: 

C/new = V (^fUoM + E ^bend,^ (10) 

where P is a projector over SU(3), and f/bcnd,i are the staples attached to 
the link in the spatial directions. Two iterations of fuzzing with / = 2.5 are 
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used and then the fuzzed hnks of length one are used. The fuzzed fermionic 
fields are defined following [17]. 

We employed two types of hadronic operator for the heavy-light mesons - 
local and fuzzed. Then for the initial state of two such mesons we have 4 
basis states. If one restricts to the operators symmetric under interchange, 
then this leaves three operators, symbolically LL, LF+FL and FF. We then 
have a 3 X 3 matrix of correlations at time t between these states at initial 
and final time. From this we use a variational approach to extract the linear 
combination of operators which maximises the ground state contribution. 

4 Results 

Our results are for quenched lattices at (3 = 5.7 and we set the scale [16] from 
(string tension) =0.44 GeV (which imphes ro = 0.53 fm) using ro/a — 2.94 
to obtain 1/a — 1.10 GeV. This scale has systematic errors of at least 10% 
coming from the differences relative to experiment of different observables 
in the quenched approximation. There will also be lattice corrections which 
should be dominantly of order since we use clover improvement. Because 
of the similarity with the lattice spacing and GeV units, we present most of 
our results in lattice units with the understanding that they can be read as 
GeV to get an estimate of the physical units. Thus we are able to measure the 
strength of the interaction out to separations of ~ 8 which will correspond 
roughly to 1.4 fm. 

For the B meson itself, needed to evaluate binding energies, we follow Ref. [16] 
and use either variational analyses or a fit to all correlations over a range 
of i-values. We find that there are substantial excited state contributions 
and that a good two-state fit is possible to our correlations from 60 gauge 
configurations for 5 <t with x^/dof = 2.4/(15 — 6) yielding = 0.876(6). 
This can be contrasted with the value of 0.875(6) obtained in Ref. [16] from 
a fit for 5 < t to a larger variational basis from 20 gauge configurations. For 
a variational study, we determine the basis from using i of 3 and 4 and then 
follow the effective mass in that basis to larger t to look for a plateau which 
we find by t values of 6 and 7 - see Table 2. This gives similar results to the 
fit approach. 
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For a study of the BB system, one approach would be to use the variational 
basis found in the B meson study for each of the two B mesons. This will 
certainly be a good approach at large R when any interaction between the 
two B mesons will be very small. We shall use this basis to give an overview 
of the relative size of different contributions to the interaction. 

A more sophisticated approach would be to make a new variational study of 
the BB system itself. The spatially-symmetric sector is described by a 3 x 3 
matrix as discussed above. We find in practice that this BB optimal basis 
gives very similar results to using the B meson basis for each B meson. 

Given that a combined fit was found to be the method of choice for the B 
meson study [16], we should also investigate fits to the BB correlations. One 
problem is that if the B is described by two states, then BB will require three 
energy eigenstatcs {BB, B'B, and B'B'). This increase in parameters makes 

the fit less stable. 

In each case, we can use a bootstrap method to study the binding energy by 
using the same subsets of gauges for the BB and B studies. 

For this study, we use on-axis separations i? = 0, 1, . . . , 5 for spatial size 12^ 
and = 0, 1, . . . , 8 for 16^. We also measured the correlation for the off-axis 
separation ol R— (±1, ±1, 0) in both cases. 

Overview of results. We first discuss our results from 12^ spatial lattices in 
terms of the ratios of contributions to the uncrossed diagram for spin average 
(C/), taking the B meson basis discussed above. For the BB correlator in 
this basis divided by the square of the B correlator in the same basis, we 
find the results given in Fig. 2. This shows that, for i? > 2, we find this 
ratio to be consistent with constant versus t. This constancy implies that 
there would be no binding energy for this correlation within the errors. The 
fact that the ratio is larger than one can be explained as a consequence of 
the fluctuation with spatial location and gauge conflguration of the B meson 
correlator and the property that (c^) > (c)^ for a fluctuating quantity c. The 
relatively large statistical error on the signal of interest, the departure of the 
ratio in Fig. 2 from 1.0, is consistent with the observation that disconnected 
quark diagrams (as this is) are noisier in lattice simulations than connected 
diagrams (such as the crossed diagram). 

The ratio of the spin-averaged BB correlation from the cross diagram to that 
from the uncrossed diagram is shown in Fig. 3. The ratio is seen to increase 
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with t and to decrease with R {R = is anomalous). This t-dependence 
imphes an interaction, and we find it to decrease in relative strength with 
increasing R. 

The uncrossed spin flip correlation (C^ averaged over x, y and z) is fairly 
small as shown in Fig. 4 and has big errors. The dominant contribution to 
the spin- flip comes from the cross diagram as illustrated in Fig. 5. In both 
cases, the spin-flip correlation is poorly determined at larger R. We shall 
discuss these contributions in terms of particle exchange later. 

We flnd that the uncrossed diagram mainly contributes to the spin average, 
while the crossed diagram contributes a comparable amount to the spin-flip 
and spin average. This is easy to understand since crossing the quarks will 
cause the spin average component (siS2 — >■ s^s^ of H — H — ) to become 
H — — h which is spin flip. We also looked at the tensor interaction 
{2Cs{z) — Cs{x) — Cs{y)) but found a small and poorly determined signal. 

In the analysis presented above, the B meson ground state has been extracted 
by using the variational basis found from a study of a single B meson. It is 
not feasible to construct a pure two meson state on a lattice in Euclidean 
time since asymptotic states cannot be constructed. Rather, one can only 
construct a state with given quantum number and then extract the energy 
eigenvalues. Nevertheless, a qualitative understanding can be obtained, as 
above, by constructing approximations to the two meson state and exploring 
their correlations. 

For example if we find that the ratio Cg/Ci ^ C+Dt, and if the combinations 
Ci + fiCg and C/ -|- /2Cs correspond to two given sets of quantum numbers 
1, 2 in the BB channel, then the mass difference which is obtained from the 
^-dependence of these correlators at large t satisfies 

Ei-E2 = j^ log ^ - h){C + Dt) ^ {h - h)D (11) 

if Cs/Cj is small. Thus a hnear dependence of the ratio is expected and 
can be related to the energy difference as shown. We do indeed see evidence 
for such linear behaviours of ratios in the figures just discussed. This linear 
dependence of the spin-fiip correlation on t can also be related theoretically 
to a meson exchange interpretation, for example, and we discuss this later. 

Energy Levels. A study of correlations between lattice operators at increasing 
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t allows an analysis of energy levels. Thus by taking the appropriate com- 
bination of crossed and spin-flip contributions, the energy of BB states with 
different quantum numbers can be studied. Because the binding energies are 
found to be quite small, we have presented the foregoing discussion to show 
the quality of our lattice data. 

We present the energies for isospin and 1 light quarks for the triplet and 
singlet spin combinations, using a 3 x 3 variational basis from t of 4 and 3 to 
obtain the optimum combination for the BB ground state. Here we use C/ -|- 
Cg correlation for the triplet states with Cg the average over orientation which 
is appropriate for an S'-wave bound state and C/ — 3(7^ for the singlet states. 
The energies evaluated on a lattice include a contribution from the self-energy 
of the static source which is unphysical. Thus only energy differences have 
a physical significance and hence we concentrate especially on the binding 
energies - the difference of BB energy from twice the B energy. In the special 
case of = 0, we show the actual lattice energy values in Table 2 to allow 
us to discuss the extrapolation to large t needed to extract the ground state. 
Other results are given in Table 3 for the case of R = 3 and in Figs. 6 to 
9. We show the results from both 12^ and 16'^ spatial lattices with the same 
parameters in order to explore finite size effects. Within errors, we do not 
see significant differences in the results between spatial sizes of L = 12 and 
16, which is not unexpected since a study of the B meson using L — 8 and 
12 found [16] agreement for the energies of the ground state mesons and a 
relatively localised Bethe-Saltpctcr wave function. 

The situation at /? = is special because the two static b quarks can be 
classified under their combined colour into either an anti-triplet or a sextet. 
The former case is just that which applies to baryons with one static quark 
and these are expected to be the lightest states. Thus the BB spectrum at 
i? = is expected to reproduce these baryonic levels. As shown in Table 2, 
we find excellent agreement with the masses of baryonic states with one 
static quark which have been obtained on the lattice previously [16]. This 
is a useful cross check of our procedures for obtaining energy levels. Thus 
we find that the A5 (with light quarks of / = and in a spin singlet) is the 
lightest state. Combining with the B meson mass then gives an estimate of 
the binding energies at i? = which will agree well with those we obtain 
here - namely around 400 MeV for the Iq, Sq = (0, 0) state. We are also able 
to explore the energies of states with R — having the opposite symmetry 
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Table 2: Effective masses for B and for BB at R=0 

Iq Sq L Eefl 

i-ratio: ref[16] 









5/4 


6/5 


7/6 




B 

1/2 


1/2 


12 


0.911(3) 


0.893(3) 


0.873(6) 


0.875(6) 


BB 














1 


1 


12 


1.620(9) 


1.516(12) 


1.558(32) 


1.514(52) 






16 


1.589(10) 


1.539(18) 


1.476(35) 










12 


1.472(18) 


1.412(29) 


1.301(63) 


1.435(37) 






16 


1.458(17) 


1.420(31) 


1.348(48) 




1 





12 


1.864(22) 


1.806(44) 


1.629(111) 








16 


1.915(19) 


1.882(50) 


1.827(130) 







1 


12 


1.911(19) 


1.852(36) 


1.722(109) 








16 


1.860(18) 


1.865(43) 


2.886(134) 





- so corresponding to the sextet of colour which is symmetric (rather than 
the anti-triplet which is antisymmetric). We find these states to he higher 
in mass than the anti-triplet states by about 0.3 in lattice units as shown in 
Table 2 and to be unbound. 

Unlike on the lattice, where for static quarks, the binding energy at i? = 
can be obtained by taking the difference of the lattice baryon mass with twice 
the fattice B mass, in the continuum, in the heavy quark limit, one would 
expect that the binding of the BB system at i? = for light quark flavour of 
/ = is given by 2{Mb — mf,) — (Ma^ — m^,) where mi, is the 6-quark mass. 

Since we flnd that the variational method gives a plateau from values of 5 
and 6, as shown in Table 2, we expect that that would be a good criterion 
to use for R> 0. 

However, for the B meson case itself, it is found that our variational method 
does not achieve a plateau value for the effective mass until a t-ratio of 7 
to 6 as shown in Table 2. At these large t-values, the BB signal is very 
noisy. Since the same operators are used for the BB case as for the B meson 
alone, it is feasible that excited state contributions are dealt with similarly 
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Table 3: Binding energies for BB at R=3a 

Ig Sg L E{BB)-2E{B) 

i-ratio: 









5/4 


6/5 


7/6 


1 


1 


12 


0.027(5) 


0.019(13) 


0.012(32) 






16 


0.024(4) 


0.023(12) 


0.005(41) 








12 


0.035(10) 


0.050(30) 


0.021(90) 






16 


-0.013(13) 


0.034(32) 


-0.000(85) 


1 





12 


-0.002(10) 


-0.077(23) 


0.040(72) 






16 


-0.030(10) 


-0.003(22) 


0.069(92) 





1 


12 


-0.029(6) 


-0.040(10) 


-0.061(38) 






16 


-0.021(7) 


-0.038(12) 


-0.019(42) 



in each case, particularly for i? > where the binding energy is found to be 
very small. Thus it makes sense to study the difference (the binding energy) 
obtained from the BB effective mass at a given t-ratio and twice the B meson 
effective mass at the same t-ratio. This is plotted in Figs. 6 to 9 from the 
ratio of correlations at t-values of 5 and 6. 

To explore the consistency of the binding energy obtained in this way, we 
show in Table 3 at i? = 3 the variational effective mass differences from 
different values. This leads us to conclude that the variational effective 
mass values for the binding energy are consistent with being constant within 
errors from t- values of 5/4, the excited state contamination being smaller 
than for the total energies. For extra safety in extracting the ground state, 
we shall use the effective mass from the values of 6/5, as stated above. 

As a cross check of this procedure, we find that the binding energy is consis- 
tent with zero within errors at large R, namely R> 5. 

As one goes to nonzero R, the level ordering found ai R — would be 
expected to be retained if the dominant dynamical configuration was that 

the two heavy quarks combine to an anti-triplet. We illustrate the binding 
energies for these states analogous to the A;, in Fig. 7 and the Sf, in Fig. 6. 
We see the level ordering to persist for the smallest values of R, the binding 
disappearing at it! 0.2 fm for the E^-hke state and at it! 0.3 fm for the 
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Afe analogue. The binding for the other pair of states is shown in Figs. 8 and 
9. Here we see that the Iq, Sq={0,l) state shows a statistically significant 
binding of 40 MeV at 0.5 fm. The situation for the Ig, Sg={l,0) state is 
less clear, since the statistical fluctuation is larger, but it is consistent with a 
similar interpretation. Note that pion exchange in the cross diagram will act 
to make the 1^, Sg = (1,0) and (0,1) states lightest at large R as we discuss 
in more detail later. 

5 Discussion 

Bound states. We find binding at small R for Ig, Sq={0,0) and (1,1) and 
binding at moderate R (circa 0.5 fm) for (1,0) and (0,1). For very heavy 
quarks, this will imply binding of the BB molecules with these quantum 
numbers and L — 0. For the physically relevant case of b quarks of around 
5 GeV, the kinetic energy will not be negligible and the binding energy of 
the BB molecular states is less clear cut. One way to estimate the kinetic 
energy for the BB case with reduced mass circa 2.5 GeV is to use analytic 
approximations to the potentials we find. For example the Iq, Sq={0,0) case 
shows a deep binding at i? = which we can approximate as a Coulomb 
potential of —0.1/R in GeV units. This will give a di-meson binding energy 
of only 10 MeV. For the other interesting case, Ig, Sg={0,l), a harmonic oscil- 
lator potential in the radial coordinate of form — 0.04[1 — (r — 3)^/4] in GeV 
units leads to a kinetic energy which completely cancels the potential energy 
minimum, leaving zero binding. This harmonic oscillator approximation lies 
above our estimate of the potential, so again we expect weak binding of the 
di-meson system. 

Because of these very small values for the di-meson binding energies, we need 
to retain corrections to the heavy quark approximation to make more definite 
predictions, since these corrections are known to be of magnitude 46 MeV for 
the B system. It will also be necessary to extrapolate our light quark mass 
from strange to the lighter u, d values to make more definite predictions 
about the binding of B mesons. This is especially necessary for light meson 
exchange contributions, which discuss subsequently. 

A model for static four-quark systems is extended and fitted to our binding 
energies in Ref. [19]. As in the static case, the results point out the inad- 
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equacy of a simple two-body potential approach for describing multi-quark 
systems. Inclusion of a multi-quark interaction term interpolating between 
strong and weak coupling regimes enables reproduction of the lattice data. 

One meson exchange. The interaction responsible for the binding energy in 
the BB system can be discussed in terms of meson exchange. One simple cri- 
terion is that BB — > BB only allows natural parity exchange (such as vector 
meson exchange) while BB* ^ B*B has an unnatural parity exchange com- 
ponent as well. Here natural means that the exchanged mesons have parity 
(— 1)*^. This can be explored by viewing the diagrams of Fig. 1 as represent- 
ing a (spatially non-local) meson creation at z — zi and then annihilation 
at Z2 = Zi + R. The quantum numbers of the mesons propagating in the 
z-direction then can be determined from the Dirac structure of the effective 
creation operator. So for Cj (BB BB), we have scalar and vector mesons 
allowed (natural parity exchanges), while for Cs{z), we have pseudoscalar 
and axial (unnatural parity exchanges), while for Cs{x) and Cs{y), both ax- 
ial and vector arc allowed. From this analysis it follows that at large R, the 
correlations at fixed t behave as exp(— Mi?) with M the mass corresponding 
to the lightest meson exchange allowed. For our lattice parameters, these 
will be the pseudoscalar meson, mass 0.529(2), for Cs{z) and vector meson, 
mass 0.815(5), for C/. 

Meson exchange contributes to the uncrossed diagram with flavour singlet 
exchange only while the crossed diagram has both flavour singlet and non- 
singlet mesons exchanged. In the quenched approximation, the flavour singlet 
and flavour non-singlet mesons are degenerate. However, in full QCD, the 
flavour singlet mass is modified by quark loop effects which are not present in 
the quenched case. These effects are responsible for the 77, 77' mass splitting, 
for example. Thus to make the cleanest comparison with meson exchange, 
it is appropriate to use the flavour non-singlet mesons (tt, p etc) which con- 
tribute only to the crossed diagram. Furthermore, our determinations of the 
contributions from the uncrossed diagram are considerably more noisy, so 
this comparison with the crossed diagram alone will be a tighter test. 

Then, as shown in Fig. 10, we see evidence for an exponential decrease of 
the interaction with increasing separation R with a mass exponent consistent 
with that expected, namely, vector for Cj and pseudoscalar for Cs{z). This 
agreement with the nature of the lightest meson exchange is a confirmation 
that the arguments given above apply at modest R values. Since the lattice 
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operator which creates the meson is not at zero momentum, we expect non- 
exponential contributions to yield the expression (l/i?) exp(— Mi?) where we 
have assumed that a sum over the ^-direction is taken (so t is large: here we 
need t > (2i?/M)^/^ which is satisfied in our case). This expression is just 
the conventional Yukawa potential. 

It is possible to go further, since lattice estimates for the B*B7r coupling are 
available [22] from a study of the axial matrix element between B and B*. 
Indeed, as well as the coupling itself, this lattice study also measures the form 
factor - the spatial distribution of the coupling - which is found to be quite 
localised. So we are be able to evaluate the magnitude of the pion exchange 
contribution using the lattice pion - so affording a direct comparison. 

Now consider the interaction potential for B*B — > B*B with B* spin polari- 
sation in the 2;-direction, which has a one pion exchange component at large 
R, 

2j^2 -MR 

where g/ f is the pion coupling to quarks [5] and we use the value determined 
from the lattice [22] oi g = 0.42(8) and where / is the pion decay constant 
(132 MeV). Because we wish to compare with our lattice results with heavier 
light quarks, we use the lattice pion mass (Ma = 0.53). 

Then to compare with our best determined quantity, the ratio of the crossed 
diagram contribution to Cs{z) to the uncrossed contribution to Cj, we assume 
that the ratio is small so that a linear t-dependence is appropriate, as indeed 
is compatible with our lattice results in Fig. 5. This implies that 

qlizl ^ t_g^e2^ 
Cf 2 47r/2 R ^ > 

and we plot this for t = 5 in Fig. 10, using the parameters discussed above. 

The agreement is excellent - better than should be expected given that t = 
5 is used and the signal is only well measured for i? < 5. In particular, 
non- leading contributions will be of order 1 / (MR) which is relatively large, 
namely 1/{MR) — 0.47 at i? = 4 for our lattice pion exchange, also note 
that some non-relativistic treatments of pion exchange have an explicit non- 
leading correction factor given by 1 + 3 /{MR). This implies that we should 
not take our estimate of the magnitude of one pion exchange as more than 
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a rough guide at the i?-values we are able to measure. Furthermore, for 
consistency, we should use the lattice determination of / for our lattice pion 
mass (which corresponds to quarks with the strange mass), hence / will be 
somewhat larger (by a factor of around fx/fn — ^■^)- What our comparison 
does show, however, is that the pion exchange contribution to the binding 
can be identified reliably for /? ~ 0.5 fm. This allows the realistic pion mass 
to be used to give predictions for the physical case with more confidence 
because of the agreement we find for pions heavier than the physical case. 

Note that the pion exchange contribution is to Cs{z) only which will con- 
tribute a large tensor interaction. Thus, much as for the case of deuterium, 
this is likely to be responsible for mixing between S and D wave components 
in the di-meson bound states. Thus implications for bound states are not 
st raight forward . 

In deuson models, the analysis of the pion exchange contribution to the po- 
tential makes meson-antimeson states in most cases significantly more bound 
than meson-meson systems. The possibility is raised in Ref. [5] that B*B* 
states bound by pion exchange may exist. In such models, however, the 
small R behaviour of the potential is not reliable. As discussed above, the 
most fruitful way to use otu^ results would be to take our non-perturbative 
measurement of the binding energy at small R and to modify our meson 
exchange component at larger R to have the lighter pion mass which is phys- 
ically relevant. 

6 Conclusions 

We study the BB system at fixed separation R using static 6-quarks. We 
present evidence for deep binding at small R with the light quark configura- 
tion similar to that in the Af, and baryons - so that the heavy quarks are 
in a colour-triplet di-quark state (and the light quarks have Ig, Sq={0,0) and 
(1,1) respectively). This binding energy is 400 - 200 MeV at R = but is 
very short-ranged. This binding is essentially a gluonic effect and is rather 
insensitive to the light quark mass, as shown by studies of the static baryons 
with varying light quark masses [16]. At larger R, around 0.5 fm, we see 
evidence for weak binding when the light quarks are in the Ig, Sq={0,l) and 
(1,0) states. This can be related to meson exchange and we find evidence of 
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an interaction in the spin- dependent quark- exchange (cross) diagram which 
is compatible with the theoretical contribution from pion exchange in our 
study. Using lighter, and hence more physical, light quark masses, this ef- 
fect will be modified in a predictable way, although further lattice study is 
needed with light quark masses below those we use (namely strange) to con- 
firm this. Corrections also need to be evaluated to the heavy quark limit for 
applications to realistic 6-quarks and we need to use smaller lattice spacings 
so reaching closer to the continuum limit, together with gauge configurations 
which have the contributions from sea quarks included. 

Our results show that it is plausible that exotic bbqq di-mesons exist as states 
stable under strong interactions. With the future lattice developments de- 
scribed above, it will be possible to give a definite answer from first principles 
in QCD whether this is so. 
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Figure 2: Results for ratio of uncrossed diagram for the spin average (corre- 
sponding to Ci) for two B mesons divided by tiie product of two B meson 
correlators, versus t. The separation is (fancy square), 1 (x), (1,1,0) 
(fancy plus), 2 (diamond), 3(octagon), 4 (square) and 5 (*) in lattice units. 
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Figure 3: Results for ratio of cross diagram to uncrossed diagram for the spin 
average (corresponding to Cj) for two B mesons versus t. The separation R 
is (fancy square), 1 (x), (1,1,0) (fancy plus), 2 (diamond), 3(octagon), 4 
(square) and 5 (*) in lattice units. 
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Figure 4: Results for ratio of spinflip to spin average (corresponding to 
Cg/Cj) for the uncrossed diagram witli two B mesons versus t. The sep- 
aration R is (fancy square), 1 (x), (1,1,0) (fancy plus), 2 (diamond) and 
3 (octagon) in lattice units. 
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Figure 5: Results for ratio of spinflip cross diagram to spin average uncrossed 
diagram (corresponding to CgjCi) for two B mesons versus t. Tlie separation 
i? is (fancy square), 1 (x), (1,1,0) (fancy plus), 2 (diamond), 3(octagon), 
4 (square) and 5 (*) in lattice units. 
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Figure 6: Results for the binding energy between two B mesons with light 
quarks in (Ig, Sq)—{1,1) at separation R in units of i?o ~ 0.5fm. The hght 
quark mass used corresponds to strange quarks. Results from variational 
method using basis from t 4:3 and effective mass in that basis from t 6:5. 
Results at different spatial lattice sizes are displaced in R for legibility. 
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Figure 7: Results for the binding energy between two B mesons with hght 
quarks in {Ig, Sq) = {0,0) at separation R in units of i?o ~ 0.5fm. The hght 
quark mass used corresponds to strange quarks. Results from variational 
method using basis from t 4:3 and effective mass in that basis from t 6:5. 
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Figure 8: Results for the binding energy between two B mesons with hght 
quarks in {Ig, Sq) = {l,0) at separation R in units of i?o ~ 0.5fm. The light 
quark mass used corresponds to strange quarks. Results from variational 
method using basis from t 4:3 and effective mass in that basis from t 6:5. 
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Figure 9: Results for the binding energy between two B mesons with hght 
quarks in {Ig, Sq) = {0,l) at separation R in units of i?o ~ 0.5fm. The hght 
quark mass used corresponds to strange quarks. Results from variational 
method using basis from t 4:3 and effective mass in that basis from t 6:5. 
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Figure 10: The ratio of the crossed- diagram contributions to the spin aver- 
aged uncrossed contribution for the BB correlation at t=5. Shown are the 
crossed diagram correlation for the spin average (BB — > BB, C/, multiplied 
by two for clarity of presentation, fancy squares) and the spin-flip (BB* 
B*B. Cs{z), octagons). The meson exchange expressions, exp{—MR)/R, are 
compared with these results for Cs{z) (using pion exchange with M = 0.529, 
continuous line) and for Cj (rho exchange with M = 0.815, dotted). Note 
that the pion exchange expression is normalised as described in the text, 
whereas the rho exchange contribution has an ad hoc normalisation. 
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